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Abstract 


Introduction. Computer modeling allows engineers to make valid design decisions by accurately assessing the thermal 
characteristics of design objects. The implementation of digital twin technology in the process of designing technical 
facilities is the current direction of scientific research and development. To do this, it is necessary to develop computer 
models whose accuracy meets the requirements for digital twins. However, the scientific literature does not widely present 
the results of research aimed at implementing digital twin technology in the design process. The general issues related to 
the use of digital twins in various industries are mainly considered. Therefore, the objective of this study was the 
development of a digital model and a comparative analysis of the accuracy of calculations of thermal characteristics of 
the design object. 

Materials and Methods. The main tool for conducting the research was the methodology proposed by the authors for 
developing a computer model of thermal characteristics for the implementation of digital twin technology. The numerical 
solution was implemented through constructing a thermal model for calculating the temperature field based on the finite 
element method in the ANSYS engineering analysis system from ANSYS, Inc. (USA). For the analytical solution, a 
computer model of thermal characteristics developed on the basis of the state-space method, implemented in the ANSYS 
Twin Builder module, was used. The state-space model was matched to the behavior of the original thermal model through 
approximating the transfer function to the stepwise response of the thermal load using the time domain vector 
approximation method. Verification of the constructed analytical model was carried out in the engineering calculation 
system MATLAB from the MathWorks company (USA). The research was carried out for a 400V machine model 
manufactured by NPO “Stankostroenie” LLC, Sterlitamak (Russia). 

Results. The developed digital model makes it possible to calculate the thermal characteristics of the design object with 
high accuracy. The results of the comparative analysis showed a high degree of correspondence between the values of 
thermal characteristics obtained using the proposed digital model and the results of numerical simulation. The maximum 
error in calculating thermal characteristics did not exceed 0.1°C. 

Discussion and Conclusion. Computer modeling that combines numerical calculation methods and a scientific approach based 
on digital twin technology, provides obtaining the result as close as possible to the results of experiments. The digital model 
proposed in the study is an effective solution, since it provides performing calculations to evaluate thermal characteristics in real 
time, which is one of the most important requirements for the implementation of digital twin technology. 
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AHHOTalna 


Beedenue. KomupiotepHoe MoyesMpoBaHue NO3BOJIAeT HHDKeHepaM MIPHHUMaTb OOOCHOBAHHBIe MpOCKTHBIe peuleHHA 3a 
cueT TOUHOM OCHKM TeMJIOBbIX XapaKTepHCTHK OOBeEKTOB MpoeKTHpoBaHHA. AKTyaJIbHbIM HallpaBJIeHHeM Hay4HBbIx 
MCCeqOBAaHHH HM pa3spaOoTOK ABJIAeTCA peaH3alHA TEXHOJIOrHH UM(PpoOBbIX ABOMHUKOB B IIpolecce MmpoeKTHpoBaHHA 
TeXHHYeCKHX OOBeKTOB. JIM 3TOrTO HeOOXOAHMO pa3pa0aTbIBATb KOMIIbIOTEPHbIC MOM, TOUHOCTb KOTOPBIX 
COOTBETCTBYeT TpeOOBaHHAM, IIPeXbABJIACMBIM K UM(PpoBbIM ABOMHUKaM. OHako B Hay4HoH JMTepaType HeocTaTOUHO 
IIMpPOKO MpeCTaBJICHbI Pe3yJIbTAaTbI HCCC NOBaHHH, HalipaBJICHHbIXx Ha peasiM3al{Hio TEXHOJOFH WHPpoBbIX JBOMHUKOB 
B Ipoljecce npoekTupoBaHua. B ocCHOBHOM paCCMaTpHBaIOTCA OOMIMe BOMPOCHI, CBA3AHHbIe C IPHMeCHeEHHeM WHpoBbIXx 
J{BOMHHMKOB B pa3JIM4HbIX OTPacsAx MpOMbULeHHOcTH. IlosToMy WesbiO WaHHOTO HCcIeqOBaHHA ABHJIaCb pa3pa0oTKa 
WMppoBwoli MOAeIM U CpaBHUTeIbHbIM aHasIM3 TOUHOCTH PpacueTOB TeIMIOBbIX XAPaKTePHCTHK OOBEKTA MIpOCKTHPOBAHHA. 
Mamepuansi u memoodvi. B xKauecTBe OCHOBHOrO MHCTpyMeHTa Id MpOBeyeHHA UCCeqOBAaHHA BBICTyMaeT 
TIpeqOKeHHad ABTOPaMU MeTOANKa pa3spa0oTKH KOMIIBIOTePHOM MOC TeEIIOBbIX XapaKTePHCTHK JIA peaM3alHu 
TeXHOJOrHH WHPpoBLix ABOMHUKOB. UuceHHoe pelieHHe peaM30BaHO IyTeM MOCTpOeCHHA TeIWIOBOM MOJeIM IA 
pacueTa TeMiepaTypHOro MOJIA Ha OCHOBe MeTOa KOHCUHBIX 9JIEMEHTOB B CHCTeMe MH2KeHEpHOroO aHasiu3a «Ansys» OT 
kommaHun «Ansys Inc» (CLUA). Aa anammTuueckoro pelieHHa IpHMeHsAeTCA pa3pabOoTaHHad Ha OCHOBe MeTOsIa 
IIpOCcTpaHcTBa COCTOAHHM KOMIIbIOTepHad MOJ[elIb TEMIOBbIX XapaKTePHCTHK, peasIM30BaHHad B Mogyse «Ansys Twin 
Builder». Moyen mpoctpaHcTBa COCTOAHHM MpHBOAMTCA B COOTBETCTBHE C TOBeCHHeEM HCXOAHOM TeIWIOBOH MOesH 
IYTeM IIpHONMKeHHA MepexaTOuHON (yHKUMH K MOWaroBOMy OTKIMKy TemIOBOM Harpy3KH C IIPHMeHeHHeM MeTosa 
BeKTOPHOH alilpOKCHMalHH BO BpeMeHHOH OOsacTH. Bepudukalua NOCTpOeHHOM aHasMTHA4eCKOM MOJeJIM BBITIOHATACh 
B CHCTeMe HH>KeHepHBIX pacyeTos «MATLAB» ot komnanuu «The MathWorks» (CHIA). UccneqoBpanna npoBoqusuch 
Ia cTaHka Moyen 400V mpou3BozcTBa upexupuatua OOO «HIIO «Crankocrpoenne» r. Crepmutamak (Poccns). 
Pexynemamoi ucciedoeanua. PazspadoTaHa udpoBad MojelIb, WO3BOJIAFOWIad C BbICOKOM TOUHOCTHIO BBITIOHMTb pacueT 
TeIVIOBbIX XapaKTepHCTHK OOBeKTa TIpoeKTHpoBaHHA. PesyIbTAaTbI CpaBHUTebHOTO aHasIH3a TIOKA3bIBalOT BbICOKYIO CTeICHb 
COOTBETCTBHA 3HAYCHHM TEIMVIOBbIX XapaKTePHCTHK, MOJIYYCHHBIX C TMOMOMIbIO TIpeWIOKCHHOM WAPpoBol Moses, pesyIbTaTaM 
YMCIeHHOrO MoyesMpoBanHa. MakcuMasibHad MOrpelwHOCTb pacueta TeMIOBbIX XapaKTepHCTHK He peBbriaer 0,1°C. 
Ooécyotcodenue u 3akiio“enue. KomibroTepHoe MOeMpoBaHue, COUeTAIOee YHCJICHHbIC METObI pacdeTa U HaydHblii 
TOAXO, OCHOBaHHbIM Ha TeXHOJIOTHH UMpoBbIX JBOMHHKOB, MO3BOJIAIOT MOYYHTb pe3syIbTaT MaKCHMaJIbHO 
TIpHOJMKeHHBIM K pe3yIbTaTaM 9KcTIepHMeHToB. IIpeqoweHHad B YCCIeqOBaHHM WMdpoBat MOJeIb ABILACTCA 
93(pdbeKTHBHEIM PelIeHHeM, MOCKOJbKY MO3BOIACT BBIMIOJIHHTb pacueThI JIA OLCHKH TeIMJIOBbIX XAPaKTEPHCTHK B PerKHMe 
peasIbHOrO BpeMeHH, 4TO ABIACTCA ONHHUM U3 BarxkHeiWUx TpedoBaHHit IpH peaM3aljMu TeXHONOrHM WH@poBBIx 
JIBOMHHKOB. 


Ksnrouespie = cJIOBAa: IM @poBolrt JIBOMHUK, CIOKHBIM TeXHH4eckHh o0bexT, KOMIIBFOTepHOe MOeCIMpOBaHHe, 
aBTOMaTH3HpOBaHHOoe NpOCKTHPOBaHHe, TEMUepaTypHoe Noe, TCMIOBbIC XapaKTCPpHCTHUKU 


Baarojxapnocrn. ABTODBI BbIparxKaroT OaroqapHOcTE peak U peveH3eHTaM 3a BHHUMAaTeCJIbBHOeC OTHOMICHHE K CTAaTbe 
HW yKa3aHHble 3aMeCuaHHA, KOTOPble NO3BOJIMJIM NHOBbICHTb Ce KadeCcTBO. 


@unancuposanne. McciezoBaHve BbIMOJIHeHO Ip (PUHAHCoBOH NoszWep»xKKe MunuctepcTBa HaykKH H BBICLUeTO OOpa30BaHHA B 
paMkax TIporpaMMBI CTpaTerH4ecKoro akayleMH4eckoro sMepcTBa «I IpHoputet-2030» (cormamrenue Ne 075—15-2023-151). 
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Introduction. Computer modeling has traditionally been an effective tool for solving thermal problems at an early 
stage of designing complex technical facilities. However, solving the problem requires an accurate assessment of the 
thermal characteristics of the design object to reduce the negative effects caused by an increase in temperature [1]. At the 
same time, one of the effective tools for preliminary assessment of thermal characteristics is simulation modeling in an 
engineering analysis system based on advanced digital solutions and developments. In [2], for example, a developed 
digital twin for determining thermal characteristics was presented by the authors Jianying Xiao and Kangoo Fan. The 
principle of operation of the twin was to simulate the thermal behavior of an object through displaying and correcting 
thermal boundary conditions. The experimental results showed a high accuracy of the model (more than 95%), which was 
essential for improving the accuracy of modeling thermal characteristics and thermal optimization. Therefore, the urgent 
direction of scientific research and development in the field of modeling is the use of artificial intelligence systems [3] 
and digital twins [4]. In [3] by Haoran Yi and others, an interactive model for correcting thermal boundary conditions 
based on a neural network was proposed to improve the accuracy of the analysis of thermal characteristics. The 
experimental results showed that the accuracy of calculating the temperature field exceeded 98%, and the accuracy of 
predicting thermal deformation was 96%, which effectively increased the simulation accuracy. In addition, in [5] by 
Kurganova N.V. and others, it was noted that digital twins were often used to improve physical prototypes of complex 
technical objects, since they not only allow for the information support for the design process, but also contribute to 
effective design decisions based on developments in the field of artificial intelligence. 

One of the characteristic features of digital twin technology is that reduced-order computer models are often used for 
simulation [6]. Therefore, the development of computer models is one of the basic conditions for the implementation of 
digital twin technology [7]. Model order reduction is an effective and mathematically understandable approach to 
overcome the time constraints of multidimensional simulation models. In [8] by Mirzoev D.A. and others, for example, 
a simple analytical model of thermal fields was proposed for the development of digital counterparts of the industrial arc 
welding process. Bordachev E.V. and Lapshin V.P. [9] presented the results of mathematical modeling of the temperature 
in the tool-product contact zone under metal turning. This approach provides obtaining an accurate assessment of thermal 
characteristics corresponding to the results of numerical experiments in real time. Schréder C. and Matthias V. [10] 
presented a reduced-order model and proposed a new model balancing procedure based on the transformation of the state 
shift. In conclusion, they presented the results of a comparative analysis and the results from literature sources obtained 
through a series of numerical experiments. 

The use of a linear and time-invariant reduced-order computer model allows for fast simulation while maintaining 
high calculation accuracy [11]. When developing a computer model, approximation [12] of the transfer function is 
performed to approximate the state space model to the step response of the initial thermal model [13]. Since the step 
response of the thermal load is derived from the base thermal model, the digital model should provide the same values of 
thermal characteristics. 

However, despite the fact that recently there has been an increase in interest in the digital twin, the scientific literature 
does not widely present the results of research related to the implementation of digital solutions in the design of technical 
facilities. Based on a systematic review of the literature and thematic analysis of publications on digital twins, one of the 
key knowledge gaps associated with the development of mathematical, software and methodological support for high- 
precision computer models in the framework of the implementation of digital twin technology has been identified. 

In this regard, the objective of the study was to develop a reduced-order computer model and analyze its feasibility as 
part of a digital model for accurate calculation of thermal characteristics of complex technical design objects. 

To achieve that, it was necessary to build a thermal model and calculate the temperature field of the design object, 
generate independent step responses of the thermal load using the developed software scenario [14], implement a digital 
model for calculating thermal characteristics, determine the error of calculations obtained using numerical and analytical 
solutions, and conduct a comparative analysis of the simulation results. 

Materials and Methods. The construction of the temperature field of the modeling object was performed for a 
homogeneous isotropic body on the basis of the equation of nonstationary thermal conductivity: 


6T/0t = a° AT +q,/c, (1) 
where 7 — temperature (°C); ¢— time (s); a= Jd/(p-e) — diffusivity coefficient (m7/s); 2 — thermal conductivity 
coefficient (W/m:°C); c — specific heat capacity (J/kg-°C); p — material density (kg/m*); A — Laplace operator; 
qv — volumetric heat dissipation power (W/m’). 
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The heat flow in the heat transfer process was assumed to be equal to the amount of heat transferred through an 
arbitrary surface area S per unit time t. It is expressed by the following equation: 


== J J ana. (2) 


where Q — heat flow (W); g, — heat-flux rate (W/m7); S — surface area (m7). 

The density of the heat flux during heat transfer was determined from the formula: 

qe =0°S-(Ts-Ta), (3) 

where a — heat transfer coefficient (W/(m°C)); Ts — surface temperature (°C); T.. — ambient temperature (°C). 

In this regard, to calculate the temperature field of the modeling object according to formula (1), heat flows (2) and (3) 
were set, which determined the amount of heat passing through the surface per unit time. 

A component (Fig. 1) of a metal-cutting machine 400V model manufactured by NPO “Stankostroenie” LLC 
(Sterlitamak, Russia) in the form of a drillhead was selected as the object of modeling. 


Fig. 1. Geometric model of the object in ANSYS Design Modeler: 1 — housing; 2 — electric motor; 
3 — spindle assembly; Q — heat flows; qv — volume heat release power; qn — heat flux densities; 
a — heat transfer coefficient; T1-T4 — temperature sensors 


Since the amount of heat released mainly by the electric motor (electromagnetic losses) and the spindle assembly 
(mechanical friction losses) was taken into account for calculating the temperature field, electric motor q,, and 
spindle assembly q,2, were taken as internal heat sources, near which the corresponding heat flows Q; and Q2 were 
set. The densities of heat fluxes qn, qn2 were assigned to surfaces located near the electric motor and spindle 
assembly, respectively. 

Convection was determined by the heat transfer coefficient a taking into account the conditions of heat transfer 
(natural convection in air). Since the machine is in contact with a gaseous medium (air), the amount of heat given by the 
heated surface to the environment per time unit f, is directly proportional to the difference in temperature between surface 
Ts and medium T., depending on the area of the heat-emitting surface S (3). 

When constructing a thermal model of an object (solid) consisting of a homogeneous material (structural steel) with 
constant thermophysical properties and the presence of internal heat sources, the following initial and boundary conditions 
were assigned: 

— initial conditions took into account the fixation of a constant temperature over the entire surface of the modeling 
object (¢ = 0: T = To = const); 

— boundary conditions of the second kind were set by the heat flows of the electric motor (Q)), the spindle assembly (Q2), 
the density of the heat flux (gn1) from the electric motor to the front wall of the housing and (q,2) to the inner surfaces; 

— boundary conditions of the third kind were set by the heat transfer coefficient (a) for the surfaces located inside the 
body of the drillhead; 

— boundary conditions of the fourth kind for the contact joints of the surfaces took into account the perfect thermal 
contact and the absence of thermal resistance: 

oT 
"én 

Heat flows (Q:, Q2), heat flux densities (gn, qnz), aS well as volume heat dissipation capacity (qui, gz), Were assigned 
in accordance with well-known recommendations for metal-cutting machines [1]. Boundary conditions, as well as heat 
flows, were set for external surfaces. Therefore, it was taken into account that the interrelationships of thermal fields were 


T| 5 =T aid +g: (4) 


+0’ 


a on 


present only between the outer surfaces. 
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Differential equation (1), together with the initial and boundary conditions of the second (2), third (3) and fourth (4) kinds, 
is a mathematical formulation of the problem. The task was solved using numerical and analytical modeling methods. 

The numerical solution was performed on the basis of the finite element method in the ANSYS engineering analysis 
system, which is being developed by ANSYS Inc. (USA) and supplied by “Modeling and Digital Twins” AO, an 
authorized ANSYS distributor in Russia. The geometric model (Fig. 1) of the object was imported into the ANSYS 
Workbench project with the addition of the Transient Thermal analysis block. In the numerical solution, the simulation 
parameters and boundary conditions of the thermal model were calibrated (Table 1) to approximate the model temperature 
values to the experimental data. 

Table | 
Boundary conditions (parameters) of the thermal model 


issipati i “as Heat t fe 
Heat dissipation capacity Heat flows Heat flux densities ee 


Parameter coefficient 
qui, Wim? qu2, Wim? O1, W O2, W gu, W/m? gu2, W/m? a, W/(m?-°C) 
Value 6,500 1,000 28 15 32 18 15 


In the ANSYS Mechanical module, the initial (initial temperature Tp = 24 °C) and boundary (Table 1) conditions were 
assigned for the developed grid model, the temperature field was constructed (Fig. 2). In this case, the thermal 
conductivity coefficient was assumed to be equal to 4 = 60.5 W/(m-:°C) and was assigned as such for structural steel. 

The thermal model of the object included two contact connections for an electric motor and a spindle cartridge with a 
drillhead body, and contained 7 thermal boundary conditions. The developed grid model consisted of 16,309 elements 
and 58,527 nodes. 

The total simulation time of 21,600 seconds (6 hours) was divided into intervals (1 hour) and steps At = 360 seconds 
(6 minutes), a total of N = 60 steps within which the parameters of the thermal model (boundary conditions) were assumed 
to be constant and independent of time. 


Model (84) 
(4, Geometry 
83 Materials 
vik Coordinate Systems 
# te) Connections 
f--/) Mesh 
=) [4 Transient Thermal (85) 
v=o Initial Temperature 
lil] Analysis Settings 
2} Simplorer 
v®, Internal Heat Generation1 
y@, Internal Heat Generation2 
v&, Heat Flow1 
Vv, Heat Flow2 
JB Heat Flux1 
JB4 Heat Flux2 
vy Convection 
= Solution (B6) 
--/ 3} Solution Information 
y@® Temperature 
=} Simplorer 
ye TempProbe 1 
y® TempProbe2 
y%@ TempProbe3 
ye TempProbe4 


a 


Fig. 2. Temperature field of the object in ANSYS Mechanical 


The calculation of the temperature field was required to generate an independent step response of the thermal load, 
containing information about temperature variation over time (thermal characteristics) for each parameter (Table 1) of 
the thermal model separately: 

y:(t)=7,,(¢),@=Lk fel, (5) 
where A — number of temperature sensors; m — number of parameters of the thermal model. 

Step response generation was performed in the ANSYS Mechanical module using the Application Customization 
Toolkit (ACT) extension, which supported the implementation of user scenarios. This made it possible to develop a 
software script in the Python programming language to automatically generate a special set of files containing independent 
step responses [14]. 
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The analytical solution was performed on the basis of the state-space method through constructing a model of thermal 
characteristics according to the formula: 
x = Ax+ Bu 
. (6) 
y=Cx+Du 
where x = Ox / Ot — derivative of state vector x over time ¢; y and u — vectors of output and input data, respectively; 
A, B, C, D— matrices of constant coefficients. 


In equation (6), vector x=(x1,%,...,.xy)" contains state variables, input —_- vector 
U = (Gris Gv25 Q1, Q2, Anis Yn2, To1, To2, ...5 Tox)’ — values of the parameters of the thermal model (boundary and initial 
conditions), output vector y =(T\, T>, ..., Tx)’ — values of thermal characteristics. 


The transfer function model is expressed by the following equation: 


y,(t)= ES, (t)u,,i=Lk 
where H — matrix complex transfer function. 
Matrix transfer function H is obtained through applying the Laplace transform to formula (6), which is expressed by 


(7) 


the following equation: 
H(s)=C(sI—A)' B+D, (8) 
where s — Laplace complex variable; / — unity diagonal matrix. 

Transfer function (8) reflects the dependence of the Laplace transform of output variable Y(s) = H(s)U(s) on the 
Laplace transform of the input variable U(s) of model (6) under zero initial conditions x(0) = xo = 0. In this case, the 
dimension of the transfer function matrix H depends on the output value & = 4 and the dimension of the input value 
m = 10, which corresponds to the dimension of the initial step response (5). Therefore, model (6) is brought into line with 
the behavior of the source thermal model through approximating its transfer function (7) to step response (5) using the 
vector approximation method [15]. 

To construct the coefficient matrices of model (6), a reverse transition is performed from transfer function (8) to the 
model in the state space. In this case, transfer function (8) takes the form of the equation, whose denominator contains a 
characteristic polynomial of degree / = 4 (the order of the system), and the numerator contains a polynomial of degree 
z=1-1: 

ra 


G(s)= a. P(s)=bo + > as O(s)= aq + >, 28 (9) 


where a and b— coefficients of polynomials Q(s) and P(s), respectively. 

The roots of polynomials Q(s) and P(s) represent the poles and zeros of transfer function (8), respectively. The method 
of indefinite multipliers is applied to equation (9) to decompose each element of matrix H into elementary fractions. 
Denoting the poles of the characteristic polynomial by p;, we obtain an equation of the following form: 


G(s)=D+ 0" A (10) 


? 
i=l 8 — p; 


P(s) 


where R; = lim (s— p;)——~ — matrix of dimension (k x m); q — number of poles. 
Pi O(s) 
The rank of matrix R; is denoted by 7;, and its decomposition into the product of two matrices with the full rank of the 
column and row, respectively, is performed: 


REC BP” rank (R=, (11) 
The matrices of model (6) are diagonal, of dimension 4””, B’’”, C®", D®™ (n = 56) and contain elements that are 


obtained directly from the coefficients of transfer function (10). 
System matrix A and control matrix B contain the following coefficients: 


fp, 0 0 .. 07 b, O .. 0 
we bis 0 
0 0 0 b 0 
A= Pi2 : B- ql (12) 
0 0 PB 0 0 by Dim 
[0 0 0 Pm| [0 0 be 
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Output matrix C and feed forward matrix D contain the following coefficients: 


l ci wie Cig 0 eee 0 0 ane 0 0. ... 1 .. O 0 
x eee eee eee > = ose ,D= eee wee eae eae eae eae : 13 
0... O C(k-1)i C(k-1)q O.. 0 Oo .. 0 .. 1 0 ( ) 
| 0 eis 0 0 ees 0 Chie Chg 0... 0 a. 0 1 


This approach makes it possible to perform the transition from a model in the form of transfer function (8) to a model 
in state space (6). The transition algorithms are also implemented in the MATLAB engineering calculation system of The 
MathWorks, Inc. (USA), in the form of special functions «ss2tf()» and «tf2ss()». 

To obtain the values of thermal characteristics, the Cauchy problem is solved for a system of ordinary differential 
equations, since in formula (6), variable x is a derivative of the vector of states of the temperature field in time ¢. The 
solution to system of equations (6) is obtained using the fourth-order Runge-Kutta method. 

Verification of the constructed analytical model was performed through conducting a series of computational 
experiments using an application program developed in the MATLAB system, which included the implementation of the 
fourth-order Runge-Kutta method (Fig. 3). Computational experiments were conducted on a personal computer (AMD 
Ryzen 5 5600U processor with Radeon Graphics 2.30 GHz, RAM 16.0 GB, system type 64-bit Windows 10 Pro version 
21H2 operating system), whose characteristics were basic for modern computing technology. 

The constructed analytical model and the application of the fourth-order Runge-Kutta method for solving a 
system of differential equations made it possible to calculate the values of thermal characteristics with high accuracy 
(Fig. 3). Due to the fact that the maximum error values, i.e., the difference between the model values of thermal 
characteristics obtained using numerical (FE-Model) and analytical (LTI-Model) solutions, did not exceed 0.72 °C 
over the entire modeling interval. 
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Fig. 3. Graphs of thermal characteristics in MATLAB system 


Linear and Time-Invariant (LTT) Reduced Order Model (ROM) was developed in the ANSYS Twin Builder module 
of the ANSYS engineering analysis system based on the object's temperature field and the step response generated in the 
ANSYS Mechanical module. 

The implementation of the digital model of thermal characteristics based on a temperature field using vector 
approximation consists in the sequential execution of seven main stages. 
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Stage 1. Importing a geometric model of an object and building a thermal model in an engineering analysis system. 

Stage 2. Calculation of the object's temperature field based on the developed thermal model. 

Stage 3. Generation of an independent step response based on the results of numerical modeling of the thermal 
characteristics of the object. 

Step 4. Application of the vector approximation algorithm to obtain the poles and zeros of the transfer function of the 
state space model. 

Stage 5. Building a model of the state space based on a known transfer function. 

Stage 6. Development of a computer model of thermal characteristics. 

Stage 7. Implementation of the digital model of the object. 

Thus, the proposed digital model containing a computer model for an accurate assessment of the thermal 
characteristics of complex technical design objects was obtained through sequential completing all the above steps. 

Research Results. The developed computer model (Thermal_SG400V_SML1) contains 6 inputs and 4 outputs (Fig. 4); 
this provides identifying the relationship between the parameters of the thermal model and the values of thermal 
characteristics. The volume heat release power (qv, g.2), heat flows (Q:, Qs) and heat flux densities (gn, qn2) were taken as 
input data of the computer model. The output data were thermal characteristics (71-74). The computer model is part of the 
digital model (Fig. 4 a) of thermal characteristics implemented in the ANSYS Twin Builder module. 

The values of the parameters of the digital model, presented in the tabular form (Fig. 4 c), are fed to the input of the 
computer model (Fig. 4 a) using the STEP components. The graphic module shows the values of thermal characteristics 
(Fig. 4 5) obtained at the output of the computer model. 


Thermal SG400V_ SML 


1°C 

34. 

34. 

29.0 

24. 

0 5,000 10,000 ~—-15,000 t, 8 
b) 
1 2 3 4 5 6 
Time, s 3,600 7,200 10,800 14,400 18,000 21,600 

qn2 18.0 18.0 18.0 18.0 18.0 18.0 
qvl 6,500.0 6,500.0 6,500.0 6,500.0 6,500.0 6,500.0 
Ql 28.0 28.0 28.0 28.0 28.0 28.0 
Q2 15.0 15.0 15.0 15.0 15.0 15.0 
Ov2 1,000.0 1,000.0 1,000.0 1,000.0 1,000.0 1,000.0 
qn\ 32.0 32.0 32.0 32.0 32.0 32.0 


Fig. 4. Implementation of digital model in Ansys Twin Builder system: a — computer model of thermal characteristics; 
b — module for graphical representation of results; c — module for monitoring input data of the computer model 


During the development of the computer model, the limits of dimension from 2 to 4 orders, values of the target error 
¢ = 5x10°, and the tolerance for the zero order ¢) = 2x10 were set. The remaining parameters were set automatically, 
since the vector approximation method was as automated as possible in comparison to other methods that were supported 
in the ANSYS Twin Builder module. 

In the process of developing a computer model, the ANSYS Twin Builder module automatically generated a special 


matrix of approximation errors M ; =|| y;(t)- y ;(¢)||; in the time domain, each element of which reflected the difference 


between the values of the thermal characteristics of step response (5) and transfer function (7) of the model in the state 
space. It is expressed by the following equation: 
0.43 4.04 1.53 4.97 2.22 1.82 
0.21 2.22 2.54 1.07 1.16 1.95 
M= x10, (14) 
2.66 1.63 0.39 0.26 3.43 1.82 


1.46 3.21 0.06 0.84 2.95 3.50 
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In this case, the maximum relative error did not exceed value ¢ = 4.97 x 10°. At the same time, all other error values 
turned out to be less than the specified limit ¢ = 5 x 10>. Zero-error value in matrix (14) meant that the input was ignored 
due to a very small contribution. 

The estimation of the accuracy of the calculation of thermal characteristics using the proposed digital model was 
carried out through comparative analysis of the results obtained using numerical and analytical solutions. The comparative 
analysis was performed according to the criterion of maximum error. Maximum error ATinax, i.¢., the difference between 
the values of thermal characteristics obtained using numerical and analytical solutions for all temperature sensors, was 
calculated at each time by the formula: 

AT, 


max 


= max AT, (15) 


jalan 
where AT; = |Zj— Tj,<| — error; T,; and T;,— temperature values of the finite element and digital models, respectively 
G= 1, m); m — number of temperature sensors. 


To assess the digital model accuracy, an error calculation was performed, whose results were presented in the form of 
a surface (Fig. 5 a) and a linear graph (Fig. 5 5) of the maximum time error. The surface (Fig. 5 a) represented the 
calculated error values for each temperature sensor (dT1, dT2, dT3, dT4) individually and at each time point over the 
entire simulation interval. 
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Fig. 5. Errors in the calculation of thermal characteristics: a — error for each temperature sensor; 
b — maximum error for all temperature sensors 
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The graph (Fig. 5 b) showed that for the selected time interval, maximum errorA7ina, did not exceed 0.1 °C. The error 
in calculating the thermal characteristics for each temperature sensor separately (Fig. 5 a) did not exceed the specified 
limits either. 

Discussion and Conclusion. The obtained results of computational experiments and the conducted comparative 
analysis confirm the efficiency of the proposed digital model, which provides calculating thermal characteristics with 
high accuracy ATina. = 0.052 °C, for further analysis and identification of the thermal model. 

The temperature field of the design object is influenced by many factors, which complicates the determination 
of the nominal values of thermal boundary conditions. To solve this problem, the study first analyzes the key 
technologies involved in the implementation of the digital twin of a complex technical object, then builds a thermal 
model and a computer reduced-order model LTI ROM of the transient temperature field. The developed computer 
model is used as part of a digital model that provides obtaining an accurate assessment of the thermal characteristics 
of the design object and thereby increases the efficiency of design procedures in the process of computer-aided 
design of complex technical facilities. 

The accuracy and efficiency of calculations of the reduced-order computer model and the source full-order thermal 
model are evaluated by comparative analysis of the simulation results. The results of computational experiments show 
that, from the point of view of calculation accuracy, computer models of reduced order and finite element models of full 
order are generally comparable in accuracy, the maximum calculation error is within the acceptable range and does not 
exceed 0.1°C. 

The results obtained during computational experiments do not contradict the results presented in the sources of 
scientific literature on similar topics and allow us to conclude that the use of the proposed digital model is effective for 
evaluating the thermal characteristics of complex technical objects in real time, which is one of the most important 
conditions for the implementation of digital twin technology. 

However, changes in the temperature field of a complex technical facility still depend on many factors. Therefore, in 
further research, it is necessary to develop computer models of thermal deformations and conduct effective optimization 


algorithms based on artificial intelligence to provide the reliability of simulation results obtained using digital twins. 
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